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ABSTRACT 

A  2? 5  x  100°  region  along  the  celestial  equator  (Stripe  82)  has  been  imaged  repeatedly  from 
1998  to  2005  by  the  Sloan  Digital  Sky  Survey  (SDSS).  A  new  catalogue  of  ~4  million  light- 
motion  curves,  together  with  over  200  derived  statistical  quantities,  for  objects  in  Stripe  82 
brighter  than  r  ~21.5  has  been  constructed  by  combining  these  data  by  Bramich  et  al.  This 
catalogue  is  at  present  the  deepest  catalogue  of  its  kind.  Extracting  ~  130  000  objects  with 
highest  signal-to-noise  ratio  proper  motions,  we  build  a  reduced  proper  motion  diagram  to 
illustrate  the  scientific  promise  of  the  catalogue.  In  this  diagram,  disc  and  halo  subdwarfs 
are  well-separated  from  the  cool  white  dwarf  sequence.  Our  sample  of  1049  cool  white 
dwarf  candidates  includes  at  least  eight  and  possibly  21  new  ultracool  H-rich  white  dwarfs 
(Jeff  <  4000  K)  and  one  new  ultracool  He-rich  white  dwarf  candidate  identified  from  their 
SDSS  optical  and  UKIDSS  infrared  photometry.  At  least  10  new  halo  white  dwarfs  are  also 
identified  from  their  kinematics. 

Key  words:  catalogues  -  stars:  atmospheres  -  stars:  evolution  -  white  dwarfs. 


1  INTRODUCTION 

The  number  of  detected  cool  white  dwarfs  has  risen  dramatically 
with  the  advent  of  the  deep  all-sky  surveys,  like  the  Sloan  Digital 
Sky  Survey  (SDSS).  Nonetheless,  certain  classes  of  white  dwarfs  - 
for  example,  ultracool  white  dwarfs  and  halo  white  dwarfs  -  remain 
intrinsically  scarce. 

The  first  ultracool  (Te ^  <  4000  K)  white  dwarfs  were  discov¬ 
ered  by  Harris  et  al.  (1999)  and  Hodgkin  et  al.  (2000).  Molecu¬ 
lar  hydrogen  in  their  atmospheres  causes  high  opacity  at  infrared 
wavelengths,  producing  a  spectral  energy  distribution  with  depleted 
infrared  flux  (e.g.  Harris  et  al.  1999;  Bergeron  &  Leggett  2002). 
Six  ultracool  white  dwarfs  have  already  been  found  in  the  SDSS  on 
the  basis  of  their  unusual  colours  and  spectral  shape  (Harris  et  al. 
2001;  Gates  et  al.  2004).  Recently,  Kilic  et  al.  (2006)  used  a  com¬ 
bination  of  SDSS  photometry  and  United  States  Naval  Observatory 
‘B’  (USNO-B)  astrometry  to  double  the  still  tiny  sample  of  these 
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interesting  objects,  which  probe  the  earliest  star  formation  in  the 
Galactic  disc. 

The  halo  white  dwarfs  are  of  astronomical  interest  as  probes  of 
the  earliest  star  formation  in  the  proto-Galaxy  and  as  tests  of  the  age 
of  the  oldest  stars.  Liebert,  Dahn  &  Monet  (1998)  first  identified  six 
candidate  halo  white  dwarfs  on  the  basis  of  high  proper  motions. 
Subsequently,  Oppenheimer  et  al.  (2001)  claimed  the  discovery  of 
38  high  proper  motion  white  dwarfs,  but  it  is  unclear  whether  they 
are  halo  or  thick  disc  members  (Reid,  Sahu  &  Hawley  2001).  Harris 
et  al.  (2006)  recently  identified  a  sample  of  ~6000  cool  white  dwarfs 
from  SDSS  Data  Release  3  (DR3)  using  reduced  proper  motions 
(RPMs),  based  on  SDSS  and  USNO-B  combined  data  (Munn  et  al. 
2004).  The  sample  included  33  objects  with  substantial  tangential 
velocity  components  (>160  km  s-1)  and  so  are  excellent  halo  white 
dwarf  candidates. 

In  this  paper,  we  also  use  SDSS  data  to  identify  new  members  of 
the  cool  white  dwarf  population.  SDSS  Stripe  82  is  a  2° 5  x  100° 
strip  along  the  celestial  equator  which  has  been  repeatedly  imaged 
between  1998  and  2005.  We  exploit  the  new  catalogue  of  almost 
four  million  light-motion  curves  in  Stripe  82  (Bramich  et  al.,  in 
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preparation).  By  extracting  the  subset  of  ~  130  000  objects  with 
high  signal-to-noise  ratio  proper  motions,  we  build  a  clean  RPM 
diagram,  from  which  we  identify  1049  cool  white  dwarfs  up  to 
a  magnitude  r  ~21.5.1  Further  diagnostic  information  is  obtained 
by  combining  the  catalogue  with  near-infrared  photometry  from 
the  UKIRT  Infrared  Digital  Sky  Survey  (UKIDSS)  Data  Release  2 
(DR2)  (Warren  et  al.  2007b).  This  enables  us  to  present  new,  faint 
samples  of  the  astronomically  important  ultracool  white  dwarfs  and 
halo  white  dwarfs. 


2  SDSS  AND  UKIDSS  DATA  ON  STRIPE  82 


SDSS  is  an  imaging  and  spectroscopic  survey  (York  et  al.  2000) 
that  has  mapped  more  than  a  quarter  of  the  sky.  Imaging  data  are 
produced  simultaneously  in  five  photometric  bands,  namely  u,  g,  r,  i 
and  z  (Fukugita  et  al.  1996;  Gunn  et  al.  1998,  2006;  Hogg  et  al. 
2001 ;  Adelman-McCarthy  et  al.  2006, 2007).  The  data  are  processed 
through  pipelines  to  measure  photometric  and  astrometric  properties 
(Lupton,  Gunn  &  Szalay  1999;  Smith  et  al.  2002;  Stoughton  et  al. 
2002;  Pier  et  al.  2003;  Ivezic  et  al.  2004;  Tucker  et  al.  2006). 

SDSS  Stripe  82  covers  a  ~250deg2  area  of  sky,  consisting  of  a 
2? 5  strip  along  the  celestial  equator  from  right  ascension  —49? 5 
to  +49?  5.  The  stripe  has  been  repeatedly  imaged  between  June 
and  December  each  year  from  1998  to  2005.  Sixty-two  of  the 
total  of  134  available  imaging  runs  were  obtained  in  2005.  This 
multi-epoch  five-filter  photometric  data  set  has  been  utilized  by 
Bramich  et  al.  (in  preparation)  to  construct  the  Light-Motion-Curve 
Catalogue  (LMCC)  and  the  Higher  Level  Catalogue  (HLC).2  The 
LMCC  contains  3700  548  light-motion  curves,  extending  to  mag¬ 
nitude  21.5  in  u,  g,  r,  i  and  to  magnitude  20.5  in  z.  A  typical  light- 
motion  curve  consists  of  ~30  epochs  over  a  baseline  of  6-7  yr. 
The  root  mean  square  (rms)  scatter  in  the  individual  position  mea¬ 
surements  is  23  mas  in  each  coordinate  for  r  ^  19.0,  increasing 
exponentially  to  60  mas  at  r  =  21.0  (Bramich  et  al.,  in  preparation). 
The  HLC  includes  235  derived  quantities,  such  as  mean  magni¬ 
tudes,  photometric  variability  parameters  and  proper  motion,  for 
each  light-motion  curve  in  the  LMCC.  As  an  illustration  of  the  cata¬ 
logue’s  potential,  Fig.  1  shows  a  measured  proper  motion  in  RA  and 
Dec.  for  SDSS  J224845. 93— 005407.0,  one  of  the  ultracool  white 
dwarf  candidates  extracted  from  the  LMCC.  The  source  is  faint  (r 
=  20.56)  and  has  a  measured  proper  motion  of  204  ±  5  mas  yr-1. 

UKIDSS  is  the  UKIRT  Infrared  Deep  Sky  Survey  (Hewett  et  al. 
2006;  Hambly  et  al.  2007 ;  Irwin  et  al.,  in  preparation;  Lawrence  et  al. 
2007)  carried  out  using  the  Wide  Field  Camera  (Casali  et  al.  2007) 
installed  on  the  United  Kingdom  Infrared  Telescope  (UKIRT).  The 
survey  is  now  well  underway  with  three  European  Southern  Obser¬ 
vatory  (ESO)-wide  releases:  the  Early  Data  Release  (EDR)  in  2006 
February  (Dye  et  al.  2006),  the  Data  Release  1  (DR1)  in  2006  July 
(Warren  et  al.  2007a)  and  the  Data  Release  2  (DR2)  in  2007  March 
(Warren  et  al.  2007b).  In  fact,  UKIDSS  is  made  up  of  five  surveys. 
The  UKIDSS  Large  Area  Survey  (LAS)  is  a  near-infrared  counter¬ 
part  of  the  SDSS  photometric  survey.  The  UKIDSS  DR2  (Warren 
et  al.  2007b)  provides  partial  coverage,  consisting  of  observations 
in  at  least  one  of  YJHK  bands  (Hewett  et  al.  2006),  of  —2/3  of  Stripe 
82  to  a  depth  of  K  ~  20.1. 


1  Magnitudes  on  the  AB  system  are  used  throughout  this  paper. 

2  The  LMCC  and  the  HLC  will  be  publicly  released  as  soon  as  the  Bramich 
et  al.  (in  preparation)  paper,  currently  in  preparation,  is  published. 


Figure  1.  Measured  proper  motion  in  RA  (track  at  bottom  left,  left-hand 
y-scale)  and  Dec.  (track  at  top  left,  right-handy-scale)  for  the  ultracool  white 
dwarf  candidate  SDSS  J224845.93-005407.0. 

3  DATA  ANALYSIS 


3.1  A  RPM  diagram 


The  RPM  is  defined  as  H  =  m  +  5  log  /x  +  5,  where  m  is  the 
apparent  magnitude  and  /x  is  the  proper  motion  in  arcsec  per  year. 
Recently,  Kilic  et  al.  (2006)  have  combined  the  SDSS  Data  Release  2 
(DR2)  and  the  USNO-B  catalogues,  using  the  resulting  distribution 
of  RPMs  as  a  tracer  of  cool  white  dwarfs.  Here,  our  HLC  for  Stripe  82 
enables  us  to  construct  a  RPM  diagram  using  only  SDSS  data.  Mag¬ 
nitudes  used  in  the  RPM  diagram  are  mean  magnitudes,  calculated 
from  all  available  single-epoch  SDSS  measurements.  These  values 
therefore  differ  from  magnitudes  quoted  in  the  publicly  available 
SDSS  data  base  which  are  based  on  a  single  photometric  measure¬ 
ment.  Differences  are  usually  small,  of  the  order  of  some  hundredths 
of  a  magnitude.  Typical  proper  motion  errors  in  right  ascension  and 
declination  are  ~2masyr-1  for  r  ~17  and  ~3masyr-1  for  r  ~21 
which  is  slightly  better  than  that  of  Kilic  et  al.  (2006)  despite  the 
much  shorter  time  baseline  used  for  calculating  proper  motions  in 
the  HLC.  Additionally,  the  shorter  time  baseline  with  many  position 
measurements  has  the  advantage  of  reducing  the  contamination  due 
to  mismatches  of  objects  with  larger  proper  motions.  The  Stripe  82 
photometric  catalogue  extends  approximately  1.5  mag  deeper  than 
the  SDSS-DR2/USNO-B  catalogue,  corresponding  to  a  factor  of  2 
in  distance.  Taking  into  account  the  different  sky  coverage  of  the 
two  catalogues,  the  resultant  volume  over  which  white  dwarfs  may 
be  detected  in  the  HLC  is  ~60  per  cent  of  that  accessible  in  the 
SDSS-DR2/USNO-B  catalogue. 

In  the  left-hand  panel  of  Fig.  2,  we  present  the  RPM  diagram  for 
all  131  398  objects  in  Stripe  82  that  meet  the  following  criteria:  (i) 
the  light  curve  consists  of  at  least  nine  epochs  in  the  g,  r,  i  and  z 
filters,  (ii)  the  object  is  classified  as  stellar  by  the  SDSS  photometric 
analysis  in  at  least  80  per  cent  of  the  epochs,  (iii)  the  proper  motion 
is  measured  with  a  \J Ax2  ^  9.  In  the  case  of  high  proper  motion 
objects  (/x  ^  50  mas  yr-1),  the  latter  requirement  is  relaxed  to  allow 
a  Ax2  ^  6  measurement.  The  delta  chi- squared  A/2  of  the 
proper  motion  fit  is  defined  as 


A/2  =  X 


2 

a,  con 


Za,lin  +  Xs  ,con  Xs,l  in’ 


(1) 


where  x2on  *s  the  chi-squared  of  the  a/8  measurements  for  a  model 
that  includes  only  a  mean  position  and  xjL  is  the  chi-squared  of  the 
a/8  measurements  for  a  model  that  includes  a  mean  position  and  a 
proper  motion.  This  A/2  statistic  follows  a  chi-square  distribution 
with  two  degrees  of  freedom.  A  relatively  high  Ax2  threshold  was 
adopted  in  order  to  keep  distinct  stellar  populations  in  the  RPM 
diagram  cleanly  separated. 
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Figure  2.  Left-hand  panel:  RPM  diagram  for  Stripe  82.  Plotted  magnitudes  are  inverse-variance  weighted  mean  light  curve  magnitudes  without  any  correction 
for  reddening.  The  sequences  corresponding  to  white  dwarfs,  Population  I  disc  dwarfs  and  Population  II  main-sequence  subdwarfs  are  well  distinguishable. 
Right-hand  panel:  zoom  on  the  white  dwarf  region  of  the  RPM  diagram.  The  thin  blue  line  separates  the  white  dwarfs  (blue  dots  on  the  left-hand  panel)  from 
the  dwarf  and  subdwarf  populations  (black  dots  on  the  right-hand  panel).  Previously  spectroscopically  confirmed  white  dwarfs  are  also  shown:  WDs  of  spectral 
type  DA  (338  in  total)  as  cyan  diamonds,  WDs  of  spectral  type  DB  (12  in  total)  as  black  diamonds,  WDs  of  spectral  type  DC  (29  in  total)  as  green  diamonds 
and  WDs  of  any  other  spectral  type  (50  in  total)  as  magenta  diamonds  (Harris  et  al.  2003;  McCook  &  Sion  2003;  Kleinman  et  al.  2004;  Carollo  et  al.  2006; 
Eisenstein  et  al.  2006;  Kilic  et  al.  2006;  Silvestri  et  al.  2006). 
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The  object  sample  was  then  matched  with  the  UKIDSS  DR2 
catalogue  using  a  search  radius  of  4.0  arcsec.  Approximately 
70  per  cent  of  the  sample  had  at  least  one  near-infrared  detection. 
The  unmatched  fraction  results  from  a  combination  of  the  incom¬ 
plete  coverage  of  Stripe  82  in  the  UKIDSS  DR2  and  a  proportion 
of  the  faintest  SDSS  objects  possessing  infrared  magnitudes  below 
the  UKIDSS  LAS  detection  limits. 

The  left-hand  panel  of  Fig.  2  shows  three  distinct  sequences  of 
stars,  namely  Population  I  disc  dwarfs,  Population  II  main-sequence 
subdwarfs  and  disc  white  dwarfs.  Given  that  our  primary  target 
population  consists  of  nearby  white  dwarfs,  the  magnitudes  and 
colours  used  throughout  the  paper  have  not  been  corrected  for  the 
effects  of  Galactic  reddening.  Adopting  a  boundary  between  the 
subdwarfs  and  white  dwarfs  defined  by  Hr  >  2.68(g  —  i)  +  15.21 
for  g  —  /  <  1.6  and  Hr  >  10.0(g  —  i )  +  3.5  for  g  —  i  >  1.6  produces 
a  sample  of  1049  candidate  white  dwarfs  (see  the  right-hand  panel 
of  Fig.  2).  446  of  the  white  dwarf  candidates  possess  at  least  one 
detection  in  the  near-infrared  UKIDSS  DR2.  The  discrimination 
boundary  is  similar  to  that  employed  by  Kilic  et  al.  (2006)  and 
lies  in  a  sparsely  populated  region  of  the  RPM  diagram,  where  the 
individual  magnitude  and  proper-motion  errors  are  small.3  As  a 
consequence,  the  population  of  white  dwarfs  is  not  too  sensitive  to 
the  precise  details  of  the  separation  boundary  adopted.  In  this  border 
region,  a  small  contamination  from  the  subdwarfs  is  still  possible 
(Kilic  et  al.  2006). 

In  the  left-hand  panel  of  Fig.  3,  we  show  the  overlap  in  colour 
space  between  the  cooler  end  of  the  white  dwarf  locus  with  the 
main-sequence  star  locus.  One  can  see  that  it  is  impossible  to  dis- 

3  Typical  error  in  Hr  is  0.25  and  in  g  —  i  is  0.02. 


tinguish  cooler  white  dwarfs  from  the  main- sequence  stars  based 
on  colour  information  only,  and  proper  motions  for  these  cool  ob¬ 
jects  are  needed.  The  effectiveness  of  our  RPM  selection  is  evident 
from  the  location  of  429  spectroscopically  confirmed  white  dwarfs 
from  a  number  of  sources  (Harris  et  al.  2003;  McCook  &  Sion 
2003;  Kleinman  et  al.  2004;  Carollo  et  al.  2006;  Eisenstein  et  al. 
2006;  Kilic  et  al.  2006;  Silvestri  et  al.  2006)  in  the  right-hand  panels 
of  Figs  2  and  3.  The  previously  confirmed  white  dwarfs  populate 
practically  only  the  upper  half  of  the  white  dwarf  sequence  in  the 
RPM  diagram  and  the  hotter  end  of  the  white  dwarf  locus  in  the 
colour-colour  diagram.  Confirmed  cooler  white  dwarfs  are  rare  and 
almost  all  originate  from  Kilic  et  al.  (2006)  where  a  similar  detec¬ 
tion  method  was  adopted.  The  HLC  is  the  faintest  existing  photo¬ 
metric/astrometric  catalogue  which  enables  us  to  trace  hotter  and 
thus  intrinsically  brighter  white  dwarfs  to  further  distances  and  also 
to  significantly  enlarge  the  number  of  known  cooler  and  thus  intrin¬ 
sically  fainter  white  dwarfs. 

Fig.  4  is  a  copy  of  the  right-hand  panel  of  Fig.  2  with  the  model  loci 
for  pure  H  atmosphere  white  dwarfs  overplotted  (Bergeron,  Leggett 
&  Ruiz  2001).  For  clarity,  only  loci  for  logg  =  8.0  and  9.0,  using 
typical  disc  and  halo  tangential  velocities  VT  =  30  and  160kms_1, 
are  shown.  In  reality,  the  log  g  values  for  white  dwarfs  range  from 
~7.0  to  ~9.5  and  the  tangential  velocities  in  the  disc  span  from  ~20 
to  ~40kms-1  which  results  in  a  strong  overlap  between  different 
model  loci.  Because  of  this  degeneracy,  it  is  impossible  to  determine 
unique  white  dwarf  physical  properties  only  from  their  positions  in 
the  RPM  diagram.  In  the  colour-colour  diagrams  of  Fig.  5,  the 
model  loci  for  pure  H  and  He  atmosphere  white  dwarfs  are  plotted 
for  a  range  of  log  g  values  (Bergeron  et  al.  2001).  White  dwarfs  for 
different  logg  values  have  very  similar  colours.  Colour  differences 
are  well  measurable  only  at  the  very  cool  end  of  the  model  tracks. 
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Figure  3.  Left-hand  panel:  colour-colour  diagrams  for  all  the  objects  plotted  in  the  left-hand  panel  of  Fig.  2.  White  dwarfs  (candidates  and  confirmed)  are 
overplotted  as  blue  dots.  The  overlap  between  the  cooler  end  of  the  white  dwarf  locus  and  the  main-sequence  star  locus  is  clearly  seen.  Right-hand  panel:  zoom 
on  the  white  dwarf  region  of  the  colour-colour  diagram.  The  main-sequence  star  locus  is  this  time  plotted  in  a  light  grey  colour.  Spectroscopically  confirmed 
white  dwarfs  are  shown,  using  the  same  colour  coding  as  in  the  right-hand  panel  of  Fig.  2.  From  this  plot,  one  can  see  that  our  new  white  dwarf  candidates  are 
typically  cooler  objects. 


This  holds  for  white  dwarfs  with  either  H  or  He  atmospheres.  It 
is  also  not  obviously  feasible  to  distinguish  between  the  H  and  He 
atmosphere  white  dwarfs.  In  some  parts  of  the  parameter  space, 
there  is  complete  degeneracy. 


3.2  x2  analysis 


In  order  to  compare  the  available  magnitude  information  for  each 
white  dwarf  candidate  with  the  best  matching  model  (Bergeron  et  al. 
2001),  a  minimal  normalised  x2  value  was  calculated: 


X 


2 


i= 1 


rrij  -m(reff,  lo gg)j  -  D 

^  err,/ 


(2) 


where  ra*  are  measured  magnitudes  for  a  given  object,  m(Tet f,  log  g)t 
its  corresponding  model  predictions  and  raenv-  measured  magnitude 
errors.4  The  best  solution  was  sought  on  a  grid  of  three  free  param¬ 
eters,  namely  effective  surface  temperature  (Te ff),  surface  gravity 
(log  g )  and  distance  modulus  ( D ).  For  each  of  the  1049  white  dwarf 
candidates,  the  four  SDSS  magnitudes,  g,  r,  i  and  z,  were  used  in 
the  fit.  Whenever  any  of  the  UKIDSS  magnitudes  were  available, 
the  fit  was  repeated  with  this  additional  information  included.  Only 
a  small  subset  of  70  white  dwarf  candidates  has  as  yet  all  three  /,  H 
and  UKIDSS  magnitudes  measured.5  Both  H  and  He  atmosphere 
models  were  compared  with  the  data.  In  the  attempt  to  better  un¬ 
derstand  the  degeneracy  effects,  the  best  and  also  the  second  best 
solutions  were  recorded. 

The  statistical  analysis  of  the  obtained  normalised  x2  results  was 
performed  first  on  the  complete  white  dwarf  sample,  using  only 
SDSS  magnitude  information.  In  all  cases,  the  normalised  x2  dis¬ 
tributions  are  strongly  asymmetric.  Typically,  fits  of  the  H  models 
perform  better  than  the  He  ones.  For  the  H  models,  the  normalised 


4  Magnitude  errors  smaller  than  0.01  were  replaced  with  this  minimal  value 
in  the  x2  calculation. 

5  Y UKIDSS  magnitude  was  not  included  in  the  x2  analysis  due  to  the  lack 
of  white  dwarf  models  for  this  waveband. 


X2  distribution  peaks  at  ~1.5  and  has  a  width  of  ~3  while  for  the 
He  models  the  distribution  peak  is  at  ~2.5  and  the  width  is  ~3.5. 
Differences  between  the  two  fits  are,  however,  usually  too  small  to 
clearly  distinguish  between  the  two,  especially  in  the  temperature 
range  where  the  H  and  He  models  are  degenerate. 

White  dwarf  atmosphere  models  observed  in  broad-band  colours 
are  also  degenerate  for  different  surface  gravity  values.  In  order  to 
estimate  how  strongly  this  degeneracy  affects  the  results  of  the  fit, 
we  compared  the  calculated  best  normalised  x2  solutions  with  the 
second  best  ones.  From  the  differences  in  the  calculated  x2  values, 
it  is  usually  not  possible  to  completely  discard  the  second  best  solu¬ 
tion.  For  instance,  in  the  H  model  case,  the  second  best  normalised 
X2  distribution  peaks  at  ~3  with  a  distribution  width  of  ~3.  Nev¬ 
ertheless,  for  ~60  per  cent  of  the  white  dwarf  candidates  both  best 
and  second  best  solution  predict  the  same  effective  temperature.  For 
the  remaining  white  dwarfs,  both  temperatures  usually  differ  from 
each  other  only  by  ^5  per  cent,  which  typically  corresponds  to  two 
adjacent  bins  on  the  discrete  grid  of  available  models.  Available  pho¬ 
tometric  information  is  thus  sufficient  to  reliably  estimate  the  white 
dwarf  effective  temperature.  Not  surprisingly,  the  fit  results  of  the 
surface  gravity  (and  as  a  consequence  also  distance)  are  much  less 
certain.  The  typical  log  g  difference  between  the  best  and  second 
best  solution  lies  in  the  range  from  0.5  to  1.5. 

We  used  the  subsample  of  70  white  dwarfs  with  complete 
UKIDSS  photometry  available  to  show  how  the  UKIDSS  data  im¬ 
prove  the  constraints  on  the  model  fits.  The  normalised  x2  of  the 
fit  using  the  complete  SDSS  and  UKIDSS  magnitude  information 
is  typically  smaller  by  0.5  than  the  value  from  the  fit  using  SDSS 
data  only.  In  85  per  cent  of  the  cases,  the  temperature  predictions 
are  the  same,  otherwise  the  difference  is  only  one  step  on  the  model 
temperature  grid,  which  confirms  the  reliability  of  the  temperature 
estimation.  Moreover,  in  85  per  cent  of  the  cases  also  the  surface 
gravity  solutions  are  the  same  and  in  the  remaining  15  per  cent  of 
the  cases  the  best  SDSS+UKIDSS  solution  often  corresponds  to  the 
second  best  SDSS  solution.  This  might  indicate  that  the  UKIDSS 
data  help  in  breaking  the  degeneracy  in  the  log  g  parameter.  When 
this  is  not  the  case,  the  normalised  x2  are  always  larger  (>  10). 
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Figure  4.  As  in  the  right-hand  panel  of  Fig.  2,  the  zoom  on  the  white  dwarf  region  of  the  RPM  diagram  is  shown.  White  dwarf  H  atmosphere  models  (Bergeron 
et  al.  2001)  for  log  g  =  8.0  (upper  of  the  dashed/dotted  lines)  and  log  g  =  9.0  (lower  of  the  dashed/dotted  lines),  using  typical  disc  and  halo  tangential  velocities 
Vj  =  30kms_1  (dashed  lines)  and  Vj  =  160kms_1  (dotted  lines)  are  overplotted.  Larger  red  circles  indicate  H-rich  ultracool  white  dwarf  candidates,  larger 
violet  circles  He-rich  ultracool  white  dwarf  candidates  and  smaller  magenta  circles  possible  H-rich  ultracool  white  dwarf  candidates  (see  detailed  explanation 
in  text).  Halo  white  dwarf  candidates  are  marked  as  light  green  larger  circles  and  possible  halo  white  dwarf  candidates  as  dark  green  smaller  circles  (see  detailed 
explanation  in  text).  Empty  circles  of  any  colour  mark  objects  that  were  already  previously  spectroscopically  observed,  full  circles  mark  new  detections. 


70  white  dwarfs  is  still  too  small  a  number  to  draw  definite  con¬ 
clusions,  however,  it  seems  that  the  inclusion  of  the  UKIDSS  data 
significantly  enhances  the  credibility  of  the  photometric  fits.  This 
will  be  very  important  in  the  near  future,  when  it  will  be  possible  to 
match  large  portions  of  the  SDSS/USNO-B  measurements  with  the 
rapidly  increasing  LAS  UKIDSS  data  base. 

4  CANDIDATES 
4.1  Ultracool  white  dwarfs 

The  signature  of  an  ultracool  surface  temperature  (Teff  <  4000  K)  is 
depleted  infrared  flux,  which  makes  the  infrared  (IR)  UKIDSS  mea¬ 


surements  crucial  for  the  recognition  of  the  ultracool  white  dwarf 
candidates.  Unfortunately,  only  a  small  subsample  of  the  1049  white 
dwarfs  have  at  least  one  UKIDSS  measurement  available  as  yet.  That 
is  why  we  first  examined  the  outcome  of  the  H  and  He  atmosphere 
model  fits  performed  on  the  complete  white  dwarf  sample  using 
SDSS  magnitude  information  only.  If  the  calculated  surface  tem¬ 
perature  was  smaller  than  4000  K  and  the  normalised  yf  value  of 
the  fit  was  <15,  the  object  was  added  to  the  ultracool  white  dwarf 
candidate  list.  In  the  subsample  of  white  dwarfs  with  at  least  one 
UKIDSS  measurement,  this  additional  piece  of  information  was 
then  used  as  a  confirmation  or  rejection  of  potential  ultracool  white 
dwarf  candidates.  Typically,  the  SDSS  and  UKIDSS  predictions 
agreed. 
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Figure  5.  Colour-colour  diagrams  for  all  the  white  dwarfs  selected  from  the  RPM  diagram.  H  and  He  atmosphere  white  dwarf  models  (Bergeron  et  al.  2001) 
for  different  surface  gravity  values  are  overplotted  in  black  and  grey,  respectively.  Larger  red  circles  indicate  H-rich  ultracool  white  dwarf  candidates,  larger 
violet  circle  He-rich  ultracool  white  dwarf  candidate  and  smaller  magenta  circles  possible  H-rich  ultracool  white  dwarf  candidates  (see  detailed  explanation 
in  text).  The  arrow  indicates  the  limiting  colour  value  for  objects  that  were  not  detected  in  J,  H  or  K,  respectively.  Halo  white  dwarf  candidates  are  marked  as 
light  green  larger  circles  and  possible  halo  white  dwarf  candidates  as  dark  green  smaller  circles  (see  detailed  explanation  in  text).  Empty  circles  of  any  colour 
mark  objects  that  were  already  previously  spectroscopically  observed,  full  circles  mark  new  detections. 


The  results  of  this  analysis,  graphically  presented  in  Figs  4  and 
5,  are  as  follows. 

(i)  Nine  ultracool  white  dwarf  candidates  with  H-rich  atmosphere 
among  which  seven  have  at  least  one  supportive  UKIDSS  measure¬ 
ment.  The  object  SDSS  J224206. 19+004822.7  has  already  been 
spectroscopically  confirmed  as  a  DC-type  white  dwarf  (no  strong 
spectral  lines  present,  consistent  with  an  H  or  He  atmosphere)  by 
Kilic  et  al.  (2006)  with  a  measured  surface  temperature  of  ~3400  K. 
The  object  SDSS  J233055. 19+002852.2  has  been  spectroscopically 
confirmed  also  as  a  DC-type  white  dwarf  by  the  same  authors.  Its 
measured  surface  temperature  is  ~4100K  which  is  just  above  the 
ultracool  limit. 

(ii)  One  ultracool  white  dwarf  candidate  based  on  SDSS  photom¬ 
etry  only  with  He-rich  atmosphere. 

(iii)  14  possible  ultracool  white  dwarf  candidates  H-rich.  For 
these  candidates,  the  fit  of  the  He  atmosphere  model  gives  in  fact  a 
slightly  better  result,  predicting  as  a  consequence  an  object  with  not 


so  extremely  low  surface  temperature.  Taking  into  account  small 
differences  between  the  H  and  He  colours  and  the  fact  that  H-rich 
white  dwarfs  are  in  fact  much  more  numerous  than  the  He-rich  ones 
makes  these  14  objects  still  serious  ultracool  white  dwarf  candi¬ 
dates.  The  object  SDSS  J2321 15.67+010223.8  has  been,  spectro¬ 
scopically  confirmed  as  a  DC-type  white  dwarf  by  Carollo  et  al. 
(2006).  This  object  is  also  the  only  one  where  the  fit  based  on  only 
SDSS  magnitudes  did  not  predict  an  ultracool  temperature  but  the 
inclusion  of  the  UKIDSS  K  magnitude  revealed  the  presence  of  the 
depleted  IR  flux.  There  are  three  additional  ultracool  candidates 
with  at  least  one  confirmed  UKIDSS  measurement. 

Based  on  the  positions  in  the  RPM  diagram  and  on  the  estimated 
tangential  velocities,  we  conclude  that  the  vast  majority  of  the  newly 
discovered  ultracool  white  dwarf  candidates  are  likely  members  of 
the  disc  population.  Measured  and  calculated  properties  of  the  24 
ultracool  white  dwarf  candidates  are  presented  in  Table  1.  Due  to 
the  relatively  large  temperature  errors,  estimated  to  be  at  least  250  K 
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Table  1.  Properties  of  the  ultracool  white  dwarf  candidates. 
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*5<t  magnitude  detection  limit  in  J,  H  or  K  band,  respectively. 


Table  2.  Properties  of  the  halo  white  dwarf  candidates. 


522 


S.  Vidrih  et  al. 


in 

N  GO 

X  Q 

00 


3- 


*  § 


33  § 


£>  Q 

o  S 


00  hHMhHMMMhHhHhHHMMHHMMMMhHhHMhHhHHhHMMMM 
HH  hh  HH  HH  HH  HH  HH  HH  HH  HH  HH  HH  HH  HH  HH  HH  HH  HH  HH  HH  HH  HH  HH  HH  HH  HH  HH  HH  HH  HH 


O  p  CO  SO  OO  CO  04  04 

so  co  Tt"  cn  so  ?4  h  Tt 


04  Tt  in 
co  in  d 


incqsqsqincqi>sqi>incqosi>cqo4 
d  P  co  co  d  in  in  04*  oo  in  co  so 


OO  O  O  O  O  O  O  O  O  O  O  r— , 

ooooggoooooooggg 

goinmoSSoinooomogSS 

go^mo4g£^^oo4^coo<g££ 


o  O  f-s 

o  o  2 
CO  so 


g  g  o  o  o  g  o 

1/~1  \  05  05  (*5  05 

^  2  o  in  in  2  o 

^  H  N  N  °  in 


(NinHONcoohOshtNcocOH^incjH 

oososoinso^oi^inascooooisoasoas 

hh(N(Nh(N(NhcOhMh(NMh(Nh 


^ONhCOHOcOONh^hOO 

O-OCCfNhOOCOOOhO 


inmoooomininoinoooininoooooooino 

(Ninoocoooin(N(NOsoooinos'xi-coinoososoO'xi-in(N(N 

in'tmmoo^f(^Hin'Osoco(N'OH^ini-hhTtHinin(N 


-  -  -  VJ  -  g^^xim^SSo^^oj^^^xioopcogin^^^^ 


(N  (N  h 


in 

ol- 

Ol- 

in 

in 

so 

oo 

in 

ol- 

in 

Ol- 

in 

CO 

Ol- 

in 

Ol- 

in 

ol- 

CO 

o 

in 

so 

so 

so 

l> 

in 

in 

-H 

-H 

-H 

-H 

-H 

-H 

-H 

-H 

-H 

-H 

-H 

-H 

-H 

-H 

-H 

-H 

-H 

-H 

-H 

-H 

-H 

-H 

-H 

-H 

-H 

-H 

-H 

-H 

-H 

-H 

in 

o 

04 

in 

o 

OS 

so 

CO 

o 

oo 

so 

OS 

in 

o 

in 

OS 

Ol- 

CO 

o 

o 

o 

CO 

so 

o 

OS 

in 

O' 

so 

o 

Os 

in 

o 

04 

CO 

o 

o 

o 

t> 

i> 

OS 

o 

in 

ol- 

so 

o 

so 

oo 

04 

r- 

04 

so 

(N 

d 

-H 

so 


o 

-H 


(N 

d 

-H 


o 

-H 


o 

-H 


o 

-H 


o 

-H 

o 


o 

-H 


o 

-H 


o 

-H 

so 


o 

-H 


o 

d 

-H 


o 

-H 

o 


oooooooooooooooooooooooooooooo 

-H-H-H-H-H-H-H-H-H-H-H-H-H-H-H-H-H-H-H-H-H-H-H-H-H-H-H-H-H-H 


oooooooooooooooooooooooooooooo 

dddddddddddddddddddddddddddddd 

-H-H-H-H-H-H-H-H-H-H-H-H-H-H-H-H-H-H-H-H-H-H-H-H-H-H-H-H-H-H 

o  in  os 


oooooooooooooooooooooooooooooo 

dddddddddddddddddddddddddddddd 

-H-H-H-H-H-H-H-H-H-H-H-H-H-H-H-H-H-H-H-H-H-H-H-H-H-H-H-H-H-H 

(N 


oooooooooooooooooooooooooooooo 

dddddddddddddddddddddddddddddd 

-H-H-H-H-H-H-H-H-H-H-H-H-H-H-H-H-H-H-H-H-H-H-H-H-H-H-H-H-H-H 


OO  04 

in  co  so 
^  co  o 
os  oo  in 
o  h  in 
o  O 


o  o 
+  + 
CO 
04  04 


04  in 
as  co 

in  as 
o  in 
O 

o  o 
+  + 
CO  Tj- 
co  in 
oo  oo 

^H  CO 

in  oo 


so  so  oo 
co  in  o> 


o  o  o 
o  o  o 
+  +  I 
Tt  so  oo 
so  o  in 
P  ot-'  d 
m  in  co 
as  o  o~ 

04  CO  CO 
O  O  O 

o  o  o 


so  -« 


OO  O  04  O  O 


'3s " 


I  +  + 

CO  00  oo 
in  oo  h 


®  OO  00  h 
in  in  04  04 
04  Tt  04  Tf 
co  in  in  04 
o  o  o  o 
o  o  o  o 


o 
o 
I 

^H  Tj-  CO 


oo 


co  Tt  so  as 
o  04 
oo  04  as  h 
co  >n  o 
O  O  O 

o  o  o  o 


I 

04 

in  S  S 
O  04  tJ-  04 
04  04  04  04 

o  o  Tt  m 
o  o  o  o 


.  .  CO 

o  as  04 

CO  04  — ' 
SO  04  CO 
O  Tt  O 

o  o  o 

O  O  | 


+  I 

h  in  oi  co 
oo  oo  so  in 
P  as  as 
04  04  CO 
04  o-  in  oo 
O  O  Ol-  Ol- 
04  04  04  04 

o  o  o  o 


in  oo 
P  04 
in  in 
r-  in 
04  in 
o  o 
o  o 
I  I 

CO  O 

OO  o 
in  i-h 
04  co 
co  in 
in  in 
04  04 

o  o 


04  Os 
CO  04 
CO  CO 

0-  so 

04  04 

o  o 
o  o 
I  I 

1—1  ot- 

sq 

CO  oo 
CO  04 

Tt-  as 


as  cq 
of-  04* 
O  04 
in  04 
O  CO 

o 
o  o 
+  + 
as  as 
cq  o 
P  oo 

'sj-  CO 

so 

co  in 


as  as 
so 
Tf  co 

04  CO 

CO 

O 

o  o 

+ 


oo 


I 

l> 

as 
oo  in 
O 

OO  oo 
CO  CO 
04  04 
04  04 


04  O 
in 

04  O 
04  so 
O  Tt 

o 
o  o 
I  + 

as  oo 
i>  as 
of-  so 
CO  04 

in  so 


©  2007  The  Authors.  Journal  compilation  ©  2007  RAS,  MNRAS  382,  515-525 


Cool  white  dwarfs  in  SDSS  Stripe  82 


m 

(N  GO 

X  P 
oo 


to! 

t4*  & 


=1  " 
3-  - 


< 


*  ^ 


.£>  oo 
3T  Q 
O  oo 


E  E  E  E 


oo  no  df 

3  CO*  cc 


llll 

^  ON  ”  C- 


H  h  OO  H 

NO  O  ON  00 
m  (N 


o  o  vo  o 
CO  NO  r-~  o 

VO  h  Ti  i- 


p  p 


-H  -H  -H  -H 


-H 


o  o  o  o 

-H  -H  -H  -H 

co  oo  m  o 
O  V0  t"'  t 
O  O  ON  © 


o  o  o  o 
d  d  d  d 
-H'.-H-Hdfi 

ON  CO  CO  CO 
VO  VO  oo 

d  d  on  d 


o  o  o  o 
d  d  d  d 
-H  -H  -H  -H 

H  cn  CO  ^ 
NO  tJ-  CO  OO 

on  d  on  d 


CO 


CO 


o  o  o  o 
d  d  d  d 
-H  -H  -H  -H 

CO  ^  (N  Tj- 

Cl  H  o 


OO  h  V0  ON 
CO  d  ON  NO 
—  CO  VO  — 
C  C  (N  C 


o  o  o  o 


I  I 


+  + 

CO  VO 


C-  l>  O  OO 


CO 

CO  oo 
CO  CO 
CO  CO 
CO  CO 


d-  vo 

CO  CO 
CO  CO 


NO 

o 

o 

CO 


<D  5  Ol  S 


M  *1 

Oo  ^ 
OX)  3o 


0, 

ON 

3 

*1 


oo  o  u:  oo 


d- 

o 

o 


3  d  3 


2  m  (u  m  3 


W  w 


.22  ^ 


ON 

*1  <>1  <>1 


Oo  Oo 


OX) 

o 

od 

sd 

d 


_  ==  O  04 

oo  ©>  _ 5 

CO  as 

rs)  ON 


OX) 

o 

cd 

,  sd 

s  d 

i  W 

'  o 
o 


NO 

o 

o 


§ 

E 


d 

o 


„  NO  NO  P' 

O  O  S 
O  O  O 

CO  CO  CO 

d  d  ^ 

3  3  £ 

fl  fl  g 

3  3  d 

C/3  t/3  S 

-.010101 
o  <D  <D  -d 


a 

fl 

3  .22  5  .22  3 


,  -Sh  pq 

W  •  ■ 


3  d 


-’i  “  a  a  a 


^5  Cuo  OjQ 

£Q  m  0  &o 

3  O  0X)  o 


o 

g-oc  fl  § 


Cuo  Cuo  Cuo 

0  'Ho  'Ho  'Ho 

t>0  O  O  O  _ 

o  ^  — 1  ^  •  -0  o 


C  d  d  d  ^ 


o 

d 


d  <D 
cd 

3  3 


1 1 
^  d 
©  tX) 

3  3 
3  & 

^  7d 

o  d 

3  g 

-o  5 

d  to 

3  ^ 

3  PQ 

^  <D 

§  ■S 


O 

1  'S 


_  o  ^  _ 

o  o  .22  no 


o  o 


o  o  o  o 


fl  ^ 
’53  Oh 

SI 

d  d 

d  d 
£  £ 
7d  T3 


CO  ON 

U  l\ 

to  to 

3  3 

*1  S 

d  d 


d  d 

£  r; 

cd  © 
<D  <D 


03  ^d  ,o 
^  ^  ^ 


u  s  ^ 

to  =3  to 
h"  |  Eh" 
3  3  3 
'$  2 


22  3  E 

’-H  On 

l\  l\ 

to  to 


>  > 


111 


<D  <D  <D  <D  <D 


d  t: 
^  >  > 
^d  3  -a 

3  oo  3 

3  3  3 


^  — 
f  I  |  ^hH 

d  5 

d  d 

^  is 

X3  d3 


3  3  3 
^  ^  ^ 

^  ^  <D  <D  « 

Oh  0^  o  Q  CL 

^  Jn  >>  b 


N  <  <  < 
P  Q  Q  Q 


N  < 
P  P 


<  <  < 

PPQ 


ii  ^ 

^  to 

<8  d" 

1 1 

d  c3 

>  ^ 
3  ^ 
O  3 

3  3 

<D 

<D  o, 

&  § 

li 


?l  l\  1 1  ?l  ^ 


3  | 
'S  s 
8  § 


h  h  h  h 

^d  03  03  03  o 

*tn  .”rH  .tn  ^ 

^  ^  ^  ^  3 


d  d  d 


^  £ 

cd  cd 


7d  cd 


a  3 
3  o 


<D  <D  <D  <D  d 

dd  rO  dd  rO 

^  ^  ^  ^  3 


D  <D  <D  <D  D 
OOOOO 


<  < 
P  P 


<  < 
P  P 


in 

m 

a 
3 

3  ^ 

O  H 


523 


©  2007  The  Authors.  Journal  compilation  ©  2007  RAS,  MNRAS  382,  515-525 


524  S.  Vidrih  et  al. 


(the  model  temperature  resolution  in  this  temperature  range),  some 
of  the  24  candidates  might  be  in  reality  just  above  the  ultracool 
white  dwarf  limit.  We  also  do  not  quote  the  fit  results  for  the  surface 
gravity  since  these  are  as  explained  above  too  uncertain.  It  is  pos¬ 
sible  that  some  of  the  new  ultracool  white  dwarf  candidates  close 
to  the  adopted  boundary  are  subdwarfs  instead.  All  these  open  is¬ 
sues  can  be  ultimately  resolved  only  with  spectroscopic  follow  up. 
However,  already  by  analyzing  the  SDSS  and  UKIDSS  astromet¬ 
ric/photometric  information,  it  is  certain  that  the  total  number  of 
known  ultracool  white  dwarfs  has  at  least  doubled. 

4.2  Halo  white  dwarfs 

Harris  et  al.  (2006)  present  a  sample  of  33  halo  white  dwarf  candi¬ 
dates  from  their  study  employing  the  SDSS  DR3  and  USNO-B  cat¬ 
alogues.  The  combination  of  their  larger  sky  coverage  and  brighter 
magnitude  limit  of  r  —  19.7  results  in  a  survey  volume  three  times 
that  of  our  Stripe  82  survey. 

Here,  we  adopt  the  same  criteria  of  VT  >  160  km  s-1  as  Harris 
et  al.  (2006)  to  select  halo  white  dwarf  candidates  using  the  results 
of  the  model  fit.  Unfortunately,  white  dwarf  atmosphere  models 
(Bergeron  et  al.  2001)  are  degenerate  for  different  values  of  the 
log  g  parameter  in  a  broad  temperature  range.  This  uncertainty  in 
the  calculation  of  the  surface  gravity  is  coupled  with  the  determi¬ 
nation  of  the  distance  to  the  object.  Solutions  with  larger  surface 
gravity  and  smaller  distances  are  strongly  coupled  with  those  with 
smaller  surface  gravity  and  greater  distance.  As  a  consequence,  de¬ 
termination  of  the  tangential  velocity  is  also  somewhat  uncertain. 

Therefore,  we  divide  our  halo  white  dwarf  candidates,  presented 
in  Table  2  and  plotted  in  Figs  4  and  5,  into  two  groups: 

(i)  white  dwarfs  where  the  best  and  also  second  best  fit  solu¬ 
tion  indicate  a  quickly  moving  object.  12  highly  probable  halo 
candidates  were  thus  found.  Four  were  already  previously  spectro¬ 
scopically  observed  (Kleinman  et  al.  2004;  Eisenstein  et  al.  2006) 
and  two  of  them,  namely  SDSS  J01 0207.24—003259.7  and  SDSS 
J2341 10. 13+003259.5,  were  already  included  on  the  Harris  et  al. 
(2006)  halo  candidate  list. 

(ii)  white  dwarfs  where  the  best-fitting  solution  gives  VT  > 
160kms_1  but  the  second  best- fitting  solution  does  not.  22  white 
dwarfs  were  thus  found  and  very  likely  not  all  of  them  are  halo  white 
dwarfs.  12  were  already  previously  spectroscopically  observed  by 
McCook  &  Sion  (2003),  Kleinman  et  al.  (2004)  and  Eisenstein  et  al. 
(2006). 

One  of  the  ultracool  white  dwarf  candidates,  SDSS 
J03 01 44. 09— 004439.5,  has  its  calculated  tangential  velocity 
just  below  the  threshold  limit  and  could  also  be  thus  a  member  of 
the  halo  population. 

5  CONCLUSIONS 

Bramich  et  al.  (in  preparation)  used  the  multi-epoch  SDSS  observa¬ 
tions  of  Stripe  82  to  build  the  deepest  (complete  to  r  =  21.5)  photo¬ 
metric/astrometric  catalogue  ever  constructed,  covering  ~250  deg2 
of  the  sky.  Here,  we  have  exploited  the  catalogue  to  find  new  candi¬ 
date  members  of  scarce  white  dwarf  sub-populations.  We  extracted 
the  stellar  objects  with  high  signal-to-noise  ratio  proper  motions  and 
built  a  RPM  diagram.  Cleanly  separated  from  the  subdwarf  popu¬ 
lations  in  the  RPM  diagram  are  1049  cool  white  dwarf  candidates. 

Amongst  these,  we  identified  at  least  seven  new  ultracool  H-rich 
white  dwarf  candidates,  in  addition  to  one  already  spectroscopically 
confirmed  by  Kilic  et  al.  (2006)  and  one  new  ultracool  He-rich  white 


dwarf  candidate.  The  effective  temperatures  Teff  of  these  eight  candi¬ 
dates  are  all  <4000  K.  We  also  identified  at  least  10  new  halo  white 
dwarf  candidates  with  a  tangential  velocity  of  VT  >  160kms_1. 
Spectroscopic  follow-up  of  the  new  discoveries  is  underway. 

In  the  quest  for  new  members  of  still  rare  white  dwarf  sub¬ 
populations,  precise  kinematic  information  and  even  more  the  com¬ 
bination  of  SDSS  and  UKIDSS  photometric  measurements  will  play 
an  important  role,  as  we  illustrate  here.  For  the  ultracool  white 
dwarfs  near  IR  UKIDSS  measurements  turn  out  to  be  essential  for 
making  reliable  identifications.  In  the  near  future,  the  UKIDSS  LAS 
will  cover  large  portions  of  the  sky  imaged  previously  by  SDSS.  This 
gives  bright  prospects  for  an  even  more  extensive  search  for  the  still 
very  scarce  ultracool  white  dwarfs. 
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